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I.  INTTODOCTION 

During  the  past  few  years,  there  have  been  many  advmce- 
ments  In  numerical  techniques  for  predicting  the  behavior  of 
fluid  flows.  For  exanq>le,  several  computer  models  have  been 
developed  by  Gosman,  Spalding  and  others  [1,2,3]  which  use 
the  mass,  momentum  and  energy  conservation  equations  reduced 
to  finite  difference,  nonlinear  algebraic  form.  The  develop¬ 
ment  of  reliable  con^uter  programs  of  this  ^pe  greatly  bene¬ 
fits  engineering  analysis  in  such  widely  varying  fields  as 
meteorology,  aerodynamics  and  gasdynamics. 

The  earlier  two-dimensional  computer  codes  were  based  on 
vorticity  (oi)  and  stream  function  (i|>)  [1,2,5].  This  form  of 
the  governing  equations  eliminates  pressure  emd  velocity  from 
immediate  consideration.  Pressure  is  normally  calculated 
only  after  a  converged  solution  is  obtained.  This  technique 
has  several  inherent  disadveuitages: 

1.  It  results  in  large  errors  in  the  predicted  pressure 
distributions  in  all  bub  quiescent  flow  regions  due  to  the 
higher  order  dependence  of  the  pressure  gradient  on  stream 
function  [6]. 

2.  It  is  usually  restricted  to  constant  density  flows 
or  to  flows  in  which  density  varies  only  with  temperature 
[3,6]. 

3.  The  boundary  conditions  are  difficult  to  specify 

[3,5] 

4.  Considerable  difficulty  has  been  experienced  in 
arriving  at  converged  solutions,  especially  for  non-unifonnly 
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spaced  grids  and  hi^  flow  rates  [2,5,6,71. 

5.  The  model  is  not  easily  extended  to  three- 

dimensional  flows  [3]. 

To  overcome  these  difficulties,  emjdiasis  has  been  placed  on 
developing  con^uter  codes  based  on  velocity  2md  pressure,  the 
primitive  variables. 

A  major  problem  with  any  new  computer  model  is  model 
validation.  Qie  difficulties  of  collecting  accurate  empirical 
data  are  multiplied  when  investigating  three-dimensional  and/ 
or  reacting  flows.  In  addition,  many  variables  within  these 
flows  are  not  readily  measurable  (turbulence  Intensities,  etc.). 

An  effort  to  utllze  elliptic  conputer  models  which  can 
handle  ttirbulent,  reacting,  variable  density  flows  at  high 
subsonic  and  sonic  velocities  has  been  underway  at  the  Naval 
Postgraduate  School  for  several  years.  TWo  specific  areas 
which  have  been  investigated  are  flows  in  a  turbojet  test  cell 
and  in  the  combustion  environment  of  a  solid  fuel  ramjet. 

It  is  important  to  have  the  capability  to  test  high  per¬ 
formance  jet  engines  throughout  their  operating  envelope 
under  conditions  which  approximate  installed  conditions.  This 
is  often  accomplished  in  blockhouse  type  installations  called 
turbojet  test  cells  (TJTC) .  Tlie  typical  test  cell  Incorporates 
an  inlet,  a  horizontal  test  section  and  a  vertical  exhaust 
stack.  The  engine  to  be  tested  is  normally  mounted  near  the 
center  of  the  cell  to  allow  the  development  of  a  nearly  uniform 
engine  inlet  velocity  profile.  The  engine  exhausts  Ir.t-  an 
augmentor  tube  which  entrains  additional  air  for  exhaust  gas 
cooling  and  dilution.  The  quantity  of  this  secondary  air  is 


crucial  to  proper  engine  testing  and  test  cell  performance. 

Testing  today's  high  power  and  high  mass  flow  engines  in 
these  installations  produces  a  myriad  of  noise  and  air  pollu¬ 
tion  problems.  Cell  modifications  must  often  be  made  to  mini¬ 
mize  these  problems.  This  fact  coupled  with  the  future  need 
for  larger,  more  expensive  test  cells  to  replace  obsolete  cells 
and  to  accommodate  new  generations  of  high  technology  engines, 
makes  the  development  of  reliable  modeling  methods  imperative. 
The  frequently  used  one-dimensional  models  are  not  adegtiate 
for  predicting  the  details  of  the  complicated  flows  within  a 
turbojet  test  cell  and,  therefore,  the  cells  often  do  not  per¬ 
form  to  their  designed  limits.  An  accurate  flow  model  would 
provide  a  needed  design  tool  which  could  help  prevent  costly 
design  errors  and  inprove  operating  efficiency. 

A  two-dimensional  computer  code  was  used  to  analyze 

the  flows  in  a  full  scale  and  a  subscale  turbojet  test  cell 
at  the  Naval  Postgraduate  School  [2,6].  Experimental  data 
from  the  subscale  test  cell  have  been  compared  with  confuta¬ 
tions  made  with  this  computer  model  [6].  A  primitive  variable 
(u-v-p)  confuter  model  could  improve  this  capability  by  ex¬ 
tending  it  to  specific  geometries  and  flow  rates  that  the 
ijf-u  model  was  incapable  of  predicting.  In  addition,  a  u-v-p 
model  would  more  readily  allow  variable  density  flows  to  be 
analyzed  and  should  more  accurately  predict  augmentor  pressure 
distribution.  An  application  of  a  primitive  variable  model 
to  jet  pumps  has  been  reported  by  Croft  and  Lilley  [4]. 


The  primary  objective  of  this  investigation  was  to  adapt 
a  primitive  variable  model  to  the  turbojet  test  cell  geometry 
and  flow  conditions  and  to  validate  that  model  with  data  from 
a  subscale  turbojet  test  cell.  Possible  utilization  of  the 
model  for  predicting  optimtim  placement  of  sound  suppressors/ 
jet  break-up  devices  within  the  augmentor  tube  was  also  inves¬ 
tigated. 


II.  MODEL  OVEHVIBW 


A.  INTRODUCTION 

The  coB^utar  model  used  In  this  etu^  was  adapted  from 
the  CHAMPION  2/B/FIX  con^uter  program  developed  by  Pun  and 
Spalding  [8].  CHAMPION  is  a  TWO-dimensional  Elliptic,  Fixed 
grid  oonputer  program  which  provides  a  solution  of  the  conser- 
vation  equations  for  recirculating  flows  in  finite  difference 
form. 

B.  ASSUMPTIONS 

The  flow  was  assumed  to  be  steady,  two-dimensional  and 
subsonic.  For  simplicity  the  value  of  specific  heat  (C^)  was 
assumed  to  be  constant  although  its  dependence  on  temperature 
and/or  composition  could  easily  be  included. 

A  modified  Jones-Launder  [8,10,11,12]  two  parameter  tur¬ 
bulence  model  was  incorporated  to  calculate  the  effective 
viscosity.  It  uses  five  empirical  constants  (Table  I)  and 
requires  that  two  additional  variables,  turbulence  kinetic 
energy  (K)  and  turbulence  dissipation  rate  (e) ,  be  evaluated. 
Effective  viscosity  was  calculated  using  the  formulas: 

“.ft  ’  “IM  *  “t  <!' 


where 


^t 


Cjj  p  xVe 


(2) 


The  turbulent  Prandtl  and  Schmidt  numbers  were  taken 
equal  to  unity  and,  therefore,  the  turbulent  Lewis  number  was 


€tL.  jfc  ^ 


^2  S  ®e,e£f 

1.43  1.92  0.09  1.0  1.3 


TABLE  1.  X-e  TURBULENCE  MODEL  EMPIRICAL  CONSTANTS 


unity  [9}.  The  laminar  Prandtl  number  was  also  taken  to  be 
unil^. 

C.  GOVEKNING  EQUATIONS 

The  conservation  equations  for  axi-syomietrical  flows 
with  no  tangential  variations  can  be  put  into  the  general 
form  [8]. 

(pu«)  +  (prv«)  -  It  H  -  11?  (rr^  1^  -  (3) 

convection  terms  diffusion  terms  source 

terms 

where  ^  stands  for  the  dependent  variable  (vi,v,k,e,h,  etc.) 
being  considered  (^  «>  1  for  the  continuity  equation) ,  is 

the  appropriate  effective  exchange  coefficient  for  turbulent 
flow  and  is  the  "source  term”  (Table  II) .  The  energy 

equation  in  terms  of  stagnation  enthalpy  has  no  source  terms 
since  the  turbulent  Prandtl  and  Schmidt  numbers  were  chosen 
as  unity  and  radiative  tremsport  was  neglected  [1,3].  Thus 
the  stagnation  enthalpy  is  given  by: 

h  »  h  +  (u^  +  v^)/2  +  K  (4) 

where  for  non-reacting  flows: 

h  -  Cp  T  (5) 

The  calculation  of  temperature  was  made  using  equations 
(4)  and  (5) .  Density  was  calculated  from  the  perfect  gas  law: 


p  »  P/RT 

(6) 

where 

R  -  R/il 

(7) 

TABLE  II.  GOVERNING  EQUATION  P. 


D.  OBSERVATION  OF  MASS 

On  each  radial  line  tha  mass  flow  rata  waa  calculatad 
using  tha  local  dansl^.  Iha  arror  in  mass  flow  (conparad 
to  tha  susnaation  of  "mass-in*  at  all  upstraam  boundarias)  was 
usad  to  uniformly  adjust  tha  axial  valoeity  ovar  tha  antira 
lina.  This  procass  ansurad  that  ovarall  continui^  was  satis- 
flad  on  tha  lina.  Tha  prassura  at  all  downstraam  locations 
was  than  adjustad  to  approximataly  corract  for  tha  momentum 
imbalance  created  by  tha  uniform  axial  valoci^  adjustment. 

A  "pressure  correction"  equation  was  then  solved  for  each  call 
on  tha  line.  Local  cell  velocity  (axial  and  radial)  and  pres¬ 
sure  were  then  adjusted  to  satisfy  cell-wise  continuity,  nie 
details  of  this  procedure  are  presented  in  reference  8. 

E.  BOUNDARY  CONDITIONS 
1.  Introduction 

Fixed  boundary  conditions  were  set  at  the  desired  or 
experimentally  determined  value  and  held  constant.  Specified 
gradient  boundary  conditions  were  handled  by  setting  the 
appropriate  convection/diffusion  coefficient  to  zero  in  tha 
finite  difference  equation  ("breaking  the  link")  and  than 
entering  the  appropriate  gradient  through  linearized  "falsa* 
source  terms  [8].  The  geometry  and  appropriate  boundary  con¬ 
ditions  for  tha  test  call  are  summarized  in  figure  1. 

2.  Inlet 

Although  not.  a  computar  program  limitation,  "plug 
flow"  was  assumed  at  tha  nozzle  exit  and  cell  inlet  planes. 

Tha  (secondary)  flow  inlet  of  tha  TJTC  was  recessed  approxi¬ 
mataly  0.3  meters  (figure  1)  to  allow  a  velocity  profile  to 
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develop  over  th«  length  of  the  engine.  Turbulence  kinetic 
energy  was  selected  to  be  uniform  with  a  value  %dilch  corres¬ 
ponded  to  the  approximate  turbulence  Intensity  of  the  Inlet 
flow. 

3.  Axis  of  Symmetry  and  Exit  Plane 

Radial  and  axial  gradients  mre  set  equal  to  zero  on 
the  center  line  and  exit  respectively.  The  radial  flow  velo¬ 
city  was  equated  to  zero. 

4.  Solid  Boundaries 

All  solid  boundaries  were  considered  adiabatic  with 
both  velocity  components  equal  to  zero  ("no  slip”  condition) . 

For  sln^llclty,  a  too  part  boundary  layer  was  used. 
The  border  between  the  laminar  sublayer  and  the  turbulent 
layer  was  taken  at  •  11.5  [8).  y^*  was  evaluated  at 

each  near  wall  node  (p) , 

y/ -  (8) 

Where,  for  y  *  >  11.5 
P  “ 


was  assumed  uniform  from  the  wall  to  the  near  wall  grid 
point.  Thus, 


s/v 


If  y^^  ^11.5  ,  the  wall  shear  stress 
using  the  formula: 


lam 


(Tw) 


(10) 

was  calculated 


10 


P  Up  Kp^/Vln(Ep«Cpl/Spl/Vyi„) 


where 


u”*"  ■  i  ln{E  Yp'*') 


Wall  shear  stress  was  evaluated  for  <11.5  from  the 

P 

formula: 


^lam 


Due  to  the  steep  gradients  of  properties  in  turi>ulent 
flows  near  solid  boundaries,  the  source  terms  for  K  and  e 
at  near  wall  nodes  were  expressed  in  terms  of  the  wall  shear 
stress  [1,8].  also  provides  the  boundary  condition  for 

the  u  and  v  equations.  In  the  following  equation  for  tur¬ 
bulence  dissipation  rate  (e)  at  a  near  wall  node  (p) ,  the 
length  scale  is  presumed  proportional  to  the  distance  from  the 
wall  (5) . 


Kp^/V2.43« 


Eqxiation  (13)  implies  that  the  wall  shear  stress  is  calculated 
assuming  a  linear  veloci^  profile  when  Yp^  <11.5.  A  near¬ 
wall  grid  point,  therefore,  can  lie  within  the  laminar  sublayer 
but  the  source  terms  for  K  and  e  imply  that 


nudi  greater  than  one  [10,11].  This  fact  pracludas  from 

being  significantly  leas  than  11.5. 

F.  SOLUTION  PROCEOUBE 

Five  variables  (u,v,K,c  and  h)  were  solved  using  equation 
(3)  in  finite  difference  form.  The  line  by  line  iterative 
parooedure  employed  aqpwind  differencing  and  under  relaxation  to 
promote  convergence  [8].  Pressure  (relative  to  a  selectable 
position  and  magnitude  within  the  grid)  was  obtained  from  the 
mass  conservation  imposed  on  each  radial  grid  line  and  on  each 
nodal  control  volume  as  discussed  above.  Effective  viscosity, 
temperature  and  density  were  also  obtained  as  described  above. 

A  more  detailed  explanation  of  this  procedure  can  be  found  in 
reference  8. 


III.  TDBBOJET  TEST  CELL  MOISB  CCMiSIDBIUkTIOIIS 


A.  INTRODUCTION 

Li gh thill  [13]  and  others  have  built  a  general  physical 
theory  of  aerodynaniic  sound  production.  He  has  also  special¬ 
ized  that  theory  to  the  stud^  of  turbulence  as  a  source  of 
sound  [14].  Experiments  have  shown  [14]  that  aerodynamic  noise 
generated  by  a  subsonic  cold  jet  in  the  absence  of  any  flue- 
tTiating  forces  between  the  fluid  and  solid  boundaries  exhibits 
the  following  characteristics: 

(1)  The  acoustic  power  output  varies  as  a  high  power 
(near  the  eighth)  of  the  jet  veloci^. 

(2)  A  very  broad  spectrum  with  a  peak  frequency  near 
U^/2d  ,  where  U  is  the  jet  velocity  and  d  the  nozzle  diai> 
meter. 

(3)  Almost  all  the  soimd  is  radiated  in  directions  making 
an  acute  angle  with  the  jet  axis. 

(a)  Higher  frequencies,  which  are  apparently  emitted 
mainly  from  the  vicinity  of  the  nozzle  exhaust  plane,  are  ra¬ 
diated  at  an  angle  of  45^,  or  slightly  less,  to  the  jet  axis. 

(b)  Lower  f reqtiencies ,  which  are  apparently  emitted 
from  the  jet  five  to  twen^  diameters  downstream  of  the  nozzle, 
are  radiated  at  much  smaller  angles,. 

Lighthill  postulated  [14]  that  the  high  frequency  sound 
emanates  from  the  heavily  sheared  mixing  region  near  the  nozzle 
exhaust  and  that  the  low  frequency  emanates  from  the  region  of 
more  nearly  Isotropic  turbulence  farther  downstream  in  the  core 
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of  the  jet. 

While  a  n^rlad  of  reeeardi  has  been  done  on  the  generation 
and  propagation  of  noise  from  a  free  jet,  little  research  has 
been  directed  toward  the  problem  of  aerod^amlc  noise  produced 
by  confined  jets.  The  propagation  of  sound  associated  with 
this  ^pe  of  flow  is  significantly  more  conplicated  due  to  the 
close  proximity  of  several  fixtures  which  absorb  or  reflect 
sound  in  varying  degrees.  In  a  turbojet  test  cell  the  engine 
is  tested  in  close  proximity  to  an  augmentor  tube  which  is 
used  for  mixing  and  cooling  the  exhaust  gases  before  they  are 
vented  to  the  atmosphere  through  a  stack.  Each  of  these  struc¬ 
tures  influences  the  sound  propagation  in  some  way. 

The  numerical  models  discussed  above  can  be  used  to  predict 
velocity,  pressure  and  turbulence  intensity  distributions  in 
confined  ttirbulent  flows  like  those  discussed  above.  It  may 
be  possible  to  use  the  computer  codes  to  predict  optimum  place¬ 
ment  of  sound  suppressors  in  geometrically  con^licated  flows 
like  those  in  a  turbojet  test  cell. 

An  attenpt  was  made  to  measure  the  sound  produced  within 
the  Naval  Postgraduate  School  subscale  turbojet  test  cell  aug¬ 
mentor  tube.  These  data  were  con^ared  to  the  flow  field  pre¬ 
dictions  obtained  with  the  computer  model. 

B.  EXPERIMENTAL  APPARATUS 

1.  Test  Cell 

Design  and  construction  of  the  NPS  subscale  turbojet 
test  cell  are  detailed  in  reference  15.  The  cell  is  a  one- 
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eighth  scale  model  of  a  NAS  Alameda  test  cell.  A  TF  41 
engine  was  scaled  to  one-el9d>t^  In  diameter  resulting  In  the 
mass  flow  being  scaled  by  one-slxty  fourth  to  maintain  flow 
velocities  equal  to  those  In  the  full  scale  cell.  The  test 
section  Is  enclosed  by  hinged  plexlglas  sides  to  allow  easy 
access  and  visual  monitoring  of  the  section  during  operation. 
These  sides  were  held  open  during  the  acoustic  measurements 
described  In  this  report.  The  augmentor  tube,  equipped  with 
a  straight  Inlet,  exits  the  cell  through  a  removable  wall. 

Its  downstream  end  is  normally  attached  to  a  deflector-plate- 
equipped  vertical  exhaust  stack.  For  this  experiment,  however, 
the  stack  was  removed,  allowing  the  fluid  exiting  the  augmentor 
tube  to  vent  to  the  atmosphere.  In  an  effort  to  reduce  sound 
reflection  Inside  the  test  cell,  a  one-quarter  Inch  plywood 
divider  with  three-qxxarter  Inch  foam  padding  was  extended  from 
the  augmentor  Inlet  to  the  plexlglas  sides  of  the  test  cell. 

Hie  conbustor  used  to  simulate  turbojet  tailpipe  and  nozzle 
conditions  was  supplied  with  compressed- air  from  an  Allis- 
Chalmers,  twelve  stage  axial  compressor.  The  tailpipe  was 
fitted  with  either  a  small  (2  inches  dleuneter)  or  large  (3.5 
Inches  diameter)  nozzle. 

2.  Ins  trument ation 

Several  methods  were  employed  In  an  attempt  to  measure 
the  aerodynamic  noise  generated  by  the  engine  exhaust  In  the 
subscale  test  cell.  The  proximity  of  various  structures  In¬ 
cluding  the  augmentor  tube  and  the  stack  produced  an  "acoustic, 
ally  dirty”  environment.  This  prevented  useful  data  from  being 
taken  external  to  the  augmentor  tube  regardless  of  the  pressure 
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transducer  orientation  to,  or  distance  from  the  jet  axis.  The 
method  chosen  to  obtain  these  measurements  is  illustrated  in 
figure  2. 

The  top  of  the  augmentor  tube  incorporated  a  series  of 
press vire  taps.  Several  of  these  were  modified  to  accept  a 
one-eighth  in^  directional  microphone.  The  microphone  was 
inserted  flush  with  the  inside  of  the  augmentor  tube.  Ini¬ 
tially  it  was  held  in  a  plexiglas  plug  with  a  set  screw 
(figure  2) .  The  transducer  signals  were  aoplified,  displayed 
and  recorded  using  a  Spectral  Dynamics  Corporation  model  50330 
spectrum  euialyzer  auid  a  HP  7035B  X-Y  recorder. 

After  initial  testing  em  attenpt  was  made  to  isolate  those 
frequencies  associated  with  the  physical  resonance  of  the  aug¬ 
mentor  tube  and  test  cell  from  those  generated  directly  from 
the  jet  mixing  process.  The  pressure  transducer  was  wrapped 
with  styrofoam  and  refitted  in  the  plexiglas  holder  without 
the  securing  set  screw.  To  determine  the  degree  to  which  the 
adjacent  building  was  contributing  to  the  recorded  data,  a 
plywood  sheet  was  positioned  between  the  cell  and  building. 

C.  DATA  ACQUISITION  AND  SEDUCTION 

Sufficient  data  were  recorded  for  each  test  condition  to 
allow  the  mass  flow  rates  for  the  nozzle,  test  cell  and  aug¬ 
mentor  tube  to  be  calculated.  In  addition,  a  plot  of  frequency 
vs  sound  intensl^  was  recorded  at  each  of  the  seven  augmentor 
ports.  On  subsequent  runs  with  the  modified  transducer  holder 
and  alternate  background,  data  were  recorded  under  similar 
conditions  at  augmentor  port  nvunber  four  only. 


16 


D.  COMPUTER  MODEL  PREDICTIONS 

At  atten^t  was  not  made  to  predict  the  sound  distribution 
within  the  augmentor  tube.  Rather  the  velocity,  pressvure  and 
turbulence  kinetic  energy  distributions  were  calculated  using 
the  computer  program  discussed  above.  It  was  desired  to  deter¬ 
mine  whether  these  distributions  could  be  used  as  a  first  ap¬ 
proximation  for  locating  the  peak  noise  intensities  within  the 
augmentor  tube. 


IV.  DISCUSSION  OF  RESULTS 

A.  INTRODUCTION  , 

The  purpose  of  this  study  was  to  utllze  a  primitive  verit¬ 
able  r  finite  difference  computer  program  to  analyze  the  flow 
within  a  turbojet  test  cell  2md  to  validate  the  model  with  data 
obtained  from  a  subscale  turbojet  test  cell  located  at  the 
Naval  Postgraduate  School.  Previous  work  [2,6,7]  had  accom¬ 
plished  this  task  for  a  conqputer  model  emd,  therefore, 

eiq>lrical  data  «uid  the  predictions  of  the  t)/-u  model  were  avail¬ 
able  for  comparison.  The  experimental  data  available  consisted 
of  augmentor  wall  pressure  distributions  and  radial  velocity 
profiles  along  the  length  of  the  augmentor  tube  for  low,  medium, 
and  high  engine  flow  rates. 

As  previously  indicated,  the  model  did  a  poor  job  of 

predicting  pressure  distributions  in  all  but  quiescent  flow  re¬ 
gions.  Additionally,  numerical  convergence  was  difficult  to  ob¬ 
tain  with  that  model  when  used  to  predict  high  velocity  flows 
where  compressibility  effects  are  significemt.  It  was  emticl- 
pated  that  a  primitive  variable  model  would  help  to  alleviate 
these  difficulties.  It  is  desiredsle  to  have  a  model  which  can 
be  used  to  predict  the  flow  field  for  full-throttle  engine  con¬ 
ditions  where  the  engine  exhaust  flow  is  choked. 

Adaptation  of  the  primitive  variable  model  to  the  subscale 
test  cell  geometry  required  the  use  of  several  approximations: 

(1)  In  modeling  axl-symmetrlc  flow,  the  engine  was  by 
necessity  positioned  at  the  axis  of  symmetry.  In  the  actual 
test  cell  the  engine  was  mounted  closer  to  the  deck  than  to 
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the  overhead.  It  would  be  esqpected/  therefore,  that  the 
velocity  distribution  in  the  secondary  flow  (the  flow  around 
the  engine)  would  be  somewhat  different  than  predicted. 

(2)  The  subscale  teat  cell  cross  section  was  rectangular 
while  the  engine  and  augmentor  tube  were  cylindrical.  The 
system  was  modeled  as  three  concentric  cylinders  with  cross- 
sectional  areas  equivalent  to  the  physical  system.  The  nozzle 
exit  area,  the  test  section  cross-sectional  area  and  the  eaqpir- 
ical  augmentation  ratios  and  mass  flow  rates  were  used  in  cal- 
ciilating  the  axial  inlet  velocities  used  in  both  computer  models. 

(3)  The  actual  engine  incorporated  a  converging  nozzle. 

The  engine  vaa  modeled  as  a  cylinder  with  a  diameter  equal  to 
the  actual  nozzle  exhaust  dl2uneter. 

(4)  In  the  model  the  augmentor  inlet  and  the  aft  test 
cell  wall  were  taken  to  be  flush.  The  actual  augmentor  is  often 
inserted  into  the  test  section  which  .forms  a  recirculation  zone 
above  the  augmentor  tube.  The  effects  on  the  augmentor  flow 
field  introduced  by  this  recirculation  region  have  been  shown 

to  be  mvnlmal  [7].  It  should  be  noted  that  the  \i>-u  model 
incorporated  this  recirculation  zone  and  a  flow  reducing  lip 
flange  on  the  augmentor  inlet.  When  comparisons  were  made  bet¬ 
ween  the  predictions  of  the  two  models,  the  effects  of  this  re¬ 
circulation  zone  and  the  augmentor  inlet  lip  flange  in  the  <|>-u 
model  were  minimized  by  reducing  their  dimensions  to  one  grid 
spacing. 

There  was  one  additional  difference  between  actual  test 
conditions  and  the  nvunerlcal  model.  In  a  test  cell  of  speci¬ 
fied  geometry  (cross-sectional  area,  inlet  and  stack  acoustic 
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and  pollution  troatnant  davlees,  augnontor  tube  langth,  dia¬ 
meter  and  inlet  geometry,  etc.)  the  mass  flow  rate  of  secon¬ 
dary  air  pumped  by  the  engine  exhaust  and  the  pressure  distri¬ 
bution  throughout  the  cell  are  unique  functions  of  the  engine 
flow  rate.  The  static  pressures  at  the  cell  inlet  and  stack 
exhaust  are  both  equal  to  the  local  atmospheric' pressure. 

Thus,  the  most  desirable  model  capability  would  be  the  predic¬ 
tion  of  the  cell  flow  field  as  a  function  only  of  engine  flow 
rate.  To  accomplish  this  with  a  two-dimensional  flnlte-dlffer 
ence  model  would  greatly  complicate  the  analysis  and  Increase 
the  required  computer  conqiutation  time.  In  this  investigation 
the  flow  field  was  predicted  only  for  the  test  cell  flow  from 
the  engine  inlet  to  the  exhaust  of  the  augmentor  tube.  Since 
inlet  boundary  conditions  are  required  for  the  solution  of  the 
elliptic  equations,  the  mass  flow  rate  of  secondary  air  was 
specified  at  the  experimentally  measured  valve.  The  predicted 
pressure  profiles  were  then  compared  to  the  measured  values. 

An  alternate  method  could  have  been  employed,  requiring  signi¬ 
ficantly  greater  computer  time:  iterate  on  the  inlet  mass 
flow  rate  until  a  specified  augmentor  exhaust  pressure  was 
attained. 

Three  of  the  flow  conditions  selected  for  model  valida¬ 
tion  corresponded  to  conditions  where  empirical  data  were 
availed)le.  Model  predictions  for  two  additional  conditions 
were  made  increasing  the  engine- to-augmen tor  inlet  spacing  to 
one  and  two  engine  diameters.  No  empirical  data  were  avail¬ 
able  for  the  two  latter  conditions  and,  therefore,  inlet  para¬ 
meters  were  simulated  using  en^irical  data  for  zero  engine- 


to-augmntor  entrance  spacing.  The  test  conditions  are 
stumnarized  in  Table  III. 

Figures  3  through  7  compare  predicted  axial  pressure 
distributions  and  rddial  velocity  profiles  obtained  with  the 

and  u-v-p  ooBq>uter  models.  In  addition,  the  available 
empirical  data  are  also  plotted  on  those  figures.  The  velocity 
profiles  from  both  con^uter  models  were  plotted  for  the  grid 
lines  closest  to  the  locations  of  the  experimental  data.  In 
the  cases  where  enplrlcal  data  were  not  available,  various 
representative  velocity  profiles  were  plotted  for  both  models. 
Experimental  pressxire  profiles  were  available  only  on  the  deck 
of  the  test  cell  and  along  the  top  of  the  augmentor  tube.  Thus, 
for  negative  positions  in  Figures  3  through  7  no  empirical 
data  were  available  for  direct  conparison  to  the  predicted  pro¬ 
files.  The  models  did  predict  a  very  small  pressure  drop  along 
the  test  cell  deck  in  agreement  with  experiment.  Predicted 
axial  pressure  profiles  are  presented  for  various  radial  posi¬ 
tions.  These  locations  are  given  as  fractions  of  augmentor 
tube  radius  (R  )  .  For  example,  the  pressure  distribution 
labeled  R  »  0.96  R  Indicates  that  the  distribution  is  along 
an  axial  grid  line  located  at  a  distance  96  percent  of  the  aug¬ 
mentor  radius  from  the  axis  of  symmetry.  Two  iff-u  model  axial 
pressure  profiles  are  depicted  for  each  test  condition.  One 
profile  (R  •  0.38  R  )  lies  in  the  quiescent  flow  region  bet- 
ween  the  engine  and  augmentor  and  was  previously  found  by 
Walters  [6]  to  be  the  only  location  which  produced  a  reasonable 
estimate  of  the  measured  profiles.  The  second  axial  pressure 
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TABLE  III.  TURBOJET  TEST  CELL  TEST  CONDITIGNS 


distribution  (R  •  0.13  R^)  was  locatad  at  lass  than  ona 
angina  radius  of  tha  oantar  lina.  Threa  u-v-p  modal  axial 
prassura  profilaa  ara  prasantad  for  aach  tast  condition.  Ona 
distribution  (R  *  0.04  R_)  is  naar  tha  axis  of  synnotry. 

A 

Anothar  profila  (R  >  0.28  R^)  runs  along  tha  top  of  tha 
angina  and  throixgh  tha  turbulant  mixing  ragion.  Tha  third 
profila  (R  «  0.96  R.)  is  closa  to  tha  augmantor  wall. 

B.  RESULTS  AND  DISCUSSION  -  PRESSURE  AND  VELOCITY  DISTRIBUTIONS 
1.  Test  Case  I  -  Low  Flow  Rata/Zaro  Spacing 

a.  Velocity  Profiles 

Both  models  predicted  virtimilly  Identical  velocity 
profiles  at  each  station  along  the  augmentor  tube  (Figure  3) . 
There  was  close  agreement  between  the  predicted  and  the  experi¬ 
mental  velocity  profiles  at  the  exit  of  the  augmentor.  The 
latter  result  was  expected  since  the  velocity  profile  had  be¬ 
come  fully  developed  near  the  augmentor  exit. 

b.  Pressure  Profiles 

(1)  <p-w  Model.  This  test  condition  had  the 
lowest  flow  rates  and  offered  the  best  chance  for  agreement 
with  experiment.  The  initial  pressure  drop  emd  the  pressure 
rise  were  underpredicted  in  the  outer  flow  region  (R  «  0.38  R  ) 
Nearer  the  axis  of  symmetry  (R  *  0.13  R.)  the  initial  pressure 
decrease  was  significantly  overpredicted ,  but  the  magnitude 
(from  minimum  to  maximum)  and  profile  of  the  pressure  rise  were 
in  good  agreement  with  the  experimental  wall  profile.  The  pre¬ 
dicted  pressure  curves  leveled  off  earlier  than  the  experimental 
data.  The  difficulty  in  obtaining  good  pressure  profiles  with 
this  model  are  evident.  Near  the  augmentor  exit  there  should 


b«  xMgligibl*  radial  prassuxa  variations.  The  solution  was 
converged  in  all  dependent  variables  (Table  V)  and  negligible 
pressure  variation  was  calculated  in  the  radial  direction  near 
the  augnsntor  exit.  However,  the  pressure  profiles  along  the 
two  axial  paths  (calculated  by  sumnlng  0  p/d  x) .  Ax  from  the 
inlet  end)  did  not  become  equal  near  the  augmentor  exit.  The 
latter  resulted  from  the  large  errors  in  the  predicted  profiles 
near  the  augmentor  inlet  (i.e.,  in  the  expression 

Px  “  Pref  *  ^ 
inlet) .  This  difficulty  with  the  model  has  been  previous 

ly  discussed  by  Walters  and  Netzer  [6]. 

(2)  u-v~p  Model.  Both  the  pressure  drop  and  rise 
within  the  augmentor  were  underpredicted.  However,  the  primi¬ 
tive  variable  model  more  accurately  predicted  the  location  of 
the  minimum  pressure  and  seemed  to  level  off  at  approximately 
the  same  augmentor  position  as  the  experimental  data.  Halters 
and  Netzer  [6]  have  previously  shown  that  the  models  predict 
that  mixing  is  nearly  complete  at  the  location  where  the  pres¬ 
sure  profile  levels  off  (and  simultaneously,  where  the  velocity 
distribution  approaches  a  fully  developed  profile) .  Experi¬ 
mental  data  confirmed  this  characteristic.  Figure  3  indicates 
that,  as  anticipated,  the  u-v-p  model  can  more  accurately  pre¬ 
dict  the  location  at  which  turbulent  mixing  is  complete.  It 
also  predicted  very  little  pressure  variation  with  radial 
augmentor  position  as  is  known  to  be  true  experimentally. 

2.  Test  Case  IIA  -  Medium  Flow  Rate/Zero  Spacing 
a.  Velocity  Profiles 

The  results  are  presented  in  Figure  4  and,  again 
both  computer  predictions  were  very  similar  auid  agreed  with 


(|^)  Ax  is  very  inaccurate  near  the  augmentor 


the  limited  experimental  data  at  the  augmentor  exit  plane, 
b.  Pressure  Profiles 

(1)  ip-u  Models  The  pressure  profile  nearest 
the  centerline  became  more  unrealistic  for  this  higher  flow 

s 

rate  condition.  The  initial  pressure  drop  was  greatly  exag¬ 
gerated;  however,  the  magnitude  of  the  pressure  rise  was  again 
in  good  agreement  with  experiment.  The  pressure  profile  in 
the  quiescent  flow  region  (R  «  0.38  R  )  did  not  agree  with 
the  experimental  ciirve.  Both  the  pressure  drop  and  rise  in 
the  augmentor  tube  were  underpredicted.  In  addition,  the 
minimum  pressure  point  was  predicted  to  occur  about  one  engine 
diameter  do%mstream  of  the  experimental  minimum. 

(2)  u-v-p  Model.  The  u-v-p  model  appears  to 
consistently  vnderpredlct  the  augmentor  pressure  drop  and 
rise.  The  slope  of  the  pressure  rise,  however,  was  in  reason¬ 
able  agreement  with  experiment.  As  discussed  above,  the  pre¬ 
dicted  pressure  profiles  were  determined  by  specifying  the 
inlet  mass  flow  rates  and  uniform  velocity  profiles  at  the 
engine  exhaust  and  cell  inlet.  Table  IV  presents  the  effects 
of  changes  in  the  specified  flow  rates  on  the  predicted  maxi¬ 
mum  pressure  rise.  Also  shown  in  Ted^le  IV  are  the  correspond¬ 
ing  values  predicted  using  the  one-dimensional  equations  for 
conservation  of  mass,  momemtum  and  energy.  The  latter  calcu¬ 
lations  were  made  assuming  perfect  gases,  copl'anar  engine  exit 
and  augmentor  inlet  planes,  no  losses,  and  plug  flow  for  both 
inlets.  The  1-0  calculations  do  not  correspond  directly  to 
the  2-D  calculations  (no  losses  and  plug  flow  of  secondary  air 
at  augmentor  inlet) .  However,  the  predicted  magnitudes  of 
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TABI£  IV 


^^nax  effacta  of  dianging  flow  rataa  should  bs  siadlar. 

Table  IV  indicates  that  the  2-D  and  1>D  model  predictions  for 
the  nominal  case  (Il-h)  were  in  good  agzeMieat.  However,  the 
experimental  Ap^^^  was  considerably  larger  than  predicted. 
Cases  II  E  and  II  F  show  the  predicted  effects  of  10%  changes 
in  the  engine  and  secondary  flow  rates,  respectively.  The 
predicted  1-D  effects  were  less  than  the  2-D  effects  but  in  the 
correct  direction  and  sequence.  It  is  seen  that  a  net  10%  error 

in  the  measured  augmentation  ratio  (m  from  ASMB  orifice, 

o 

m^  from  augmentor  exit  pitot  rake,  and  m^  «  m^  -  m^)  could 
produce  the  differences  between  the  experimental  and  2-D  model 
predictions.  For  exang>le,  an  larger  by  3%  and  a  total 

augmentor  flow  rate  smaller  by  3%  (which  are  within  the  experi¬ 
mental  accurary)  would  produce  a  decrease  in  augmentation  ratio 
of  approximately  8%.  Thus,  the  differences  between  the  pre¬ 
dicted  and  measured  pressure  profiles  may  not  be  as  large  as 
depicted  in  Figure  4. 

3.  Test  Caise  IIB  -  Medium  Flow  Rate/One  Engine 
Diameter  Spacing 

a.  Velocity  Profiles 

Experimental  data  were  not  available  for  this  test 
condition.  The  predicted  profiles,  axe  presented  in  Figure  5 
and  show  that  the  two  cong>uter  models  predicted  velocity  pro¬ 
files  which  were  in  close  agreement.  The  primitive  variable 
model,  however,  predicted  a  slightly  greater  initial  jet 
spreading  at  the  augmentor  entrance  and  required  a  slightly 


longer  duct  length  to  obtain  a  fully  developed  profile. 


b.  Pressurv  Profiles 

(1)  t-oi  Hodel  .  As  for  Case  IIA,  the  centemost 
pressure  distribution  greatly  exaggerated  the  initial  pressure 
drop.  In  this  case,  the  pressure  rise  also  appeared  to  be 
much  too  rapid  in  coo^ariaon  to  the  R  *  0.38  R  profile 

S 

and  the  u-v^p  profiles.  The  profile  in  the  quiescent  flow 
region  was  in  reasonable  agreement  with  the  wall  profile  ob> 
talned  with  the  u-v-p  model.  The  slopes  of  the  pressure 
rise  and  the  minimum  pressure  predicted  by  both  models  were 
nearly  Identical.  The  ^~a>  model,  however,  again  predicted 
the  minimum  pressure  to  occur  somewhat  farther  downstream  than 
did  the  u-v-p  model. 

(2)  u-v-p  Model.  For  this  case,  the  pressure 
profile  closest  to  the  augmentor  wall  had  a  significantly 
lower  minimum  pressure  them  the  other  profiles. 

4.  Test  Case  IIC  »  Medium  Flow  Rate/Two  Engine 
Diameter  Spacing 

a.  Velocity  Profiles 

No  eiiq)lrlcal  data  were  available  for  this  test 
condition.  The  computer  predictions  are  presented  in  Figure  6. 
As  for  Case  IIB,  the  primitive  variable  model  predicted  greater 
initial  jet  spreading.  In  this  case,  however,  the 
velocity  profiles  became  flat  and  the  pressure  profiles  leveled 
off  considerably  upstream  of  the  u-v-p  predicted  profiles. 

b.  Pressure  Profiles 

(1)  Model.  The  centermost  pressure  profile 

was  completely  unrealistic.  The  quiescent  region  profile  in¬ 
dicated  a  larger  pressure  drop  and  rise  and  a  minimum  pressure 
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point  farther  downstream  than  the  u-vp  model.  In  addition, 
the  model  profile  leveled  off  much  earlier. 

(2)  u-v-p  Model.  Again  the  minimum  pressure  was 
obtained  for  the  profile  closest  to  the  augmentor  wall  and 
was  located  just  inside  the  augmentor  entrance.  All  the  prlmio 
tive  variable  curves  indicate  a  more  gentle  augmentor  pressure 
rise,  leveling  about  midway  down  the  augmentor  tube. 

5.  Test  Case  III  -  High  Flow  Rate/Zero  Spacing 
a.  Velocity  Profiles 

Substantially  more  experimental  data  were  avall- 
2d>le  for  this  test  condition.  In  the  experiment  the  nozzle 
was  operated  with  a  pressure  ratio  ^Pats/^T^  less  them  cri¬ 
tical  for  one-dimensional  Isentropic  flow.  However,  using 
the  experimental  nozzle  flow  rate  and  approximating  the  nozzle 
flow  as  one-dimensional  emd  isentropic  resulted  in  a  nozzle 
exit  Mach  number  of  approximately  0.95.  This  condition  was 
imposed  on  the  models.  As  has  been  observed  for  the  previous 
test  cases,  the  velocity  profiles  for  both  models  were  quite 
similar  (Figure  7)  and  in  reasonably  good  agreement  with 
experiment.  However,  the  predicted  profiles  for  both  models 
tended  to  flatten  a  little  too  rapidly.  The  primitive  vari¬ 
able  model  again  predicted  slightly  more  initial  jet  spread¬ 
ing  at  the  augmentor  inlet  and  less  mixing  down  the  augmentor 
t\ibe  them  did  the  ip-u  model.  The  experimental  data  more 
nearly  agreed  with  the  model  at  the  augmentor  inlet  and 

with  the  u-v-p  model  downstream. 


29 


p 


Tf 


b.  Pressure  Profiles 

(1)  Model*  Again  the  center  pressure 

profile  was  completely  unrealistic.  The  quiescent  region  pres¬ 
sure  profile  underpredicted  the  augmentor  pressure  drop  and 
rise.  The  pressure  position  was  substantially  dis¬ 

placed  down  the  augmentor  tube  and  the  slope  of  the  pressure 
rise  did  not  agree  with  experiment. 

(2)  ,  u-v-p  ttodel^  All  three  u-v-p  pressure 
profiles  were  quite  similar.  They  underpredicted  both  the 
augmentor  pressure  drop  zmd  rise.  The  slope  of  the  pressure 
rise  was  in  good  agreement  with  ejqterimental  data  within  the 
first  half  of  the  augmentor  tube.  The  predicted  pressure  pro¬ 
files  from  both  the  t-w  and  u-v-p  models  leveled  off  before 
the  e:q>erimental  curves. 

C.  COMPUTATIONAL  ACCURACY  AND  REOUIRED  COMPUTER  TIME 

The  utility  of  any  computer  program  using  numerical 
methods  is  reflected  by  the  amount  of  CPU  time  required  wd 
the  ease  of  arriving  at  converged  solutions.  TeUsle  V  compares 
the  percentage  chemge  in  variable  magnitude  on  successive 
iterations  and  the  required  CPU  time. 

A  considereible  savings  in  CPU  time  was  obtained  using 
the  line-by-line  iterative  procedure  of  the  primitive  variable 
model  in  lieu  of  the  point-by-point  (Gauss-Seidel)  method 
employed  in  the  i|>-u  model. 

At  low  flow  rates,  the  convergence  was  quite  similar  for 
both  models.  However,  at  higher  flow  rates  the  u-v-p  model 
had  better  convergence  in  less  time. 
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TABLE  V.  COMPARISONS  OF  CPU  TIME  AND  (A^/(f)  *100  FOR  LAST 
CALCULATION 


"T 


/ 1  I 
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D.  COMPUTER  RELATED  PROBLEMS 

Both  models  required  that  the  proper  relaxation  para¬ 
meters  be  selected  in  order  to  obtain  convergence.  The  lade 
of  any  procedure  for  selecting  the  proper  relaxation  vi^ues 
makes  this  process  quite  tim  consuming.  Previous  research 
using  the  i(f-u  model  facilitated  the  selection  of  these 
pareuDsters  for  that  model.  It  was  found  that  the  u-v-p  model 
was  qTiite  sensitive  to  the  calculated  "pressure  corrections”. 
Obtaining  the  correct  underrelaxation  value  for  pressure 
proved  to  be  the  key  in  arriving  at  a  converged  solution  for 
the  primitive  variadcle  model. 

The  line-by-line  iterative  procedure  used  in  the  u-v-p 
model  was,  as  expected,  quite  good  in  propagating  distur¬ 
bances  downstreeua  when  sweeping  from  left  to  right.  A  down¬ 
stream  disturbance  is  propagated  upstream  by  successive  sweeps 
through  the  entire  field.  This  fact,  at  least  for  a  geometry 
incorporating  a  sudden  contraction,  made  the  convergence  de¬ 
pendent  on  the  number  of  traverses  on  each  radial  line.  An 
excessive  number  of  traverses  would  cause  divergence.  The 
nunber  of  traverses  on  each  line  was  controlled  in  two  ways. 
After  each  traverse,  residual  factors  were  calculated  for 
each  variable  and  the  largest  residual  factor  was  coiqaared  to 
a  pre-set  value.  If  the  largest  residual  was  less  than  the 
pre-set  value,  the  program  advanced  to  the  next  radial  line. 
The  program  would  also  advance  when  a  pre-set  maximum  number 
of  traverses  had  been  completed  on  any  radial  line.  To  aid 
convergence,  only  a  few  traverses  on  each  line  were  allowed 
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until  several  sweeps  through  the  field  had  allowed  the  pre¬ 
sence  of  the  contraction  wall  to  be  "felt”  upstream.  The 
number  of  traverses  on  a  line  was  then  increased.  It  was 
additionally  found  that  when  working  with  coaxial  flows  with 
radically  different  inlet  velocities,  the  normalizing  factors 
(vdiich  were  based  on  average  inlet  conditions  and  used  to 
calculate  the  residual  values  on  each  line)  resulted  in  ex¬ 
cessive  traverses  being  made  in  regions  of  high  flow  velocity. 
Repeated  calculations  on  radial  lines  which  had  already  con¬ 
verged  often  caused  divergence.  Adjusting  the  normalizing  fac¬ 
tors  downstream  of  the  engine  exit  alleviated  this  problem. 

The  primitive  variable  model  demonstrated  some  conver¬ 
gence  difficvilty  in  the  recirculation  region  adjacent  to  the 
sudden  contraction.  This  pTOblem  could  have  resulted  from 
the  relatively  large  normalizing  factors  used  in  this  local 
region  of  low  velocity  or,  as  suggested  by  Launder  and 
Spalding  [12],  it  could  have  possibly  been  due  to  the  in¬ 
adequacy  of  the  empirical  constants  (Table  I)  in  the  K-e 
model  for  the  flow  in  this  quiescent  zone. 

As  with  any  finite  difference  numerical  solution,  grid 
spacing  was  found  to  be  critical.  To  aid  convergence,  the 
grid  spacing  was  decreased  in  all  regions  where  large  pro¬ 
perty  gradients  were  expected.  Gosman,  et  al  [1]  recommended 
that  for  the  <('-u  model,  successive  spacing  should  not  in¬ 
crease  by  more  than  about  a  factor  of  1.5.  This  restriction 
was  also  employed  for  the  primitive  variable  model.  A  30  by 
30  grid  system  was  utilized  for  the  primitive  variable  model. 
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whereas  a  43  by  40  grid  was  required  for  the  ip-u  model. 

E.  SUMMhRY  OF  RESULTS  -  PRESSURE  AND  VELOCITY  DISTRIBUTIONS 


In  most  cases  the  predicted  velocity  profiles  for  both 
the  and  u-v^p  computer  models  agreed  with  each  other  and 

with  the  available  experimental  data.  The  u-v-p  model  predic¬ 
ted  more  initial  jet  spreading  at  the  augmentor  inlet,  espe¬ 
cially  in  those  cases  where  the  engine  exit  plane  was  not  flush 
with  augmentor  inlet.  Both  models  predicted  that  a  flat  velo¬ 
city  profile  was  obtained  where  the  pressure  profile  leveled 
off.  For  the  one  case  in  which  experimental  data  were  avail¬ 
able  over  the  entire  legnth  of  the  augmentor  tube,  the  pre¬ 
dicted  velocity  profiles  seemed  to  flatten  slightly  faster 
th2m  the  experimental  data. 

In  general,  the  model  demonstrated  poor  pressure 

prediction  capability.  For  low  flow  rates  the  pressure  rise 
on  the  wall  from  minimum  to  maximum  was  acciirately  predicted 
for  a  pressure  profile  calculated  near  the  axis  of  symmetry. 

For  higher  flow  rates  the  centermost  predictions  became  un¬ 
realistic  although  converged  solutions  were  obtained  for  all 
primary  variables.  For  the  predicted  pressure  profiles  along 
axial  lines  that  were  in  the  quiescent  flow  region,  the 
model  characteristically  unpredlcted  the  pressure  drop  and 
rise  along  the  augmentor  wall.  The  minimum  predicted  pres¬ 
sure  location  was  typically  somewhat  downstream  of  its  experi¬ 
mentally  determined  position. 

The  primitive  variable  model  predicted  little  radial  pres¬ 
sure  variations  except  near  the  engine  exit  plane.  It  seemed 
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to  accurately  locate  the  otinimun  pressure  position  and  to  pre¬ 
dict  the  rapid  pressure  drop  at  the  augaentor  entrance.  It 
consistently  xinderpredicted  the  pressure  drop  and  rise  in  the 
augnentor  but  tended  to  accurately  predict  the  slopes  of  the 
rising  pressure  profiles.  Eai^vex,  as  discussed  above,  reason¬ 
ably  expected  experimental  errors  could  account  for  these  dif¬ 
ferences.  Application  of  the  K-e  turbulence  model  with 
fixed  parameters  (C^,  C2,  Cj^)  to  the  test  cell  flow 

conditions  may  also  be  a  reason  for  the  lack  of  quantitative 
accuracy  in  the  predicted  pressure  profiles.  Further  e^qperi- 
mental  data  (for  example,  turbxilence  intensity  measurements) 
are  needed  to  further  validate  the  model. 

The  primitive  variable  model  appears  to  reasonably  pre¬ 
dict  the  flow  field  up  to  sonic  engine  exhaust  conditions. 

These  predictions  Include  realistic  pressure  distributions 
and  require  substantially  less  cos^uter  time  and  fewer  grid 
lines  than  the  i|;-w  model. 

F.  RESULTS  AND  DISCUSSION  -  NOISE  CONSIDERATIONS 

Computer  model  predictions  for  the  velocity,  turbulence 
kinetic  energy  and  pressure  distributions  are  presented  in 
Figures  8,  9  and  10.  The  predictions  were  for  the  two 
inches  diameter  nozzle  at  high  mass  flow  rate  (Run  1,  Table 
VI) .  The  four  curves  on  each  figure  are  for  different  radial 
positions;  jet  centerline,  just  inside  the  jet  exhaust  radius; 
just  outside  the  jet  exhaust  radius,  and  near  the  augmentor 
wall. 
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Noise  intensities  for  selected  frequency  rsnges  es  e 
function  of  augmentor  tube  position  are  displayed  in  Figures 
11,  12  and  13. 

Figures  14  through  16  and  Table  VII  compare  the  basic 
test  results  to  data  obtained  at  one  transducer  port  after 
the  microphone  was  insulated  and  the  "building  wall"  location 
was  changed. 

Figure  11  presents  the  data  for  the  small  nozzle  with  a 
high  subsonic  exhaust  velocity  (M^  <*»  0.75).  The  highest 
frequencies  were  greater  than  peak  values  expected  from  ori¬ 
fice  data  (V/2d)  but  were  apparently  generated  near  the 
nozzle  exit  as  expected.  If  the  high  frequencies  generated 
near  the  nozzle  exit  were  radiated  at  45°  to  the  axis,  the 
corresponding  peak  intensities  would  occur  between  two  and 
five  jet  diameters  downstream.  The  lowest  frequencies  peaked 
in  intensity  at  approximately  20  jet  diameters  doimstream  in 
general  agreement  with  data  from  free  jets  emitted  from  orifice 
plates  [141.  The  intensity  increased  as  the  frequency  de¬ 
creased. 

\ 

Figure  12  presents  the  data  for  the  small  nozzle  with  a 
.  very  low  exit  velocity  (M^  %  0.25).  The  data  had  similar 
characteristics  to  that  for  the  high  exit  velocity:  (1)  the 
intensity  increased  with  decreasing  frequency,  (2)  the  peak 
frequencies  were  greater  than  V/2d  and  (3)  the  lowest  fre¬ 
quencies  peaked  in  Intensity  at  approximately  20  jet  diameters 
downstream.  In  this  case,  however,  the  highest  frequencies 
peaked  in  intensity  further  downstream. 
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Pigur*  13  prasenta  data  for  the  largar  noszla  with  low 
exit  veloci^.  The  data  were  similar  to  Figure  12  except 
that  the  intensity  did  not  vary  as  consistently  with  frequency 
range.  The  peak  intensity  did  occur  for  the  lowest  frequency 
range  and  occurred  at  approximately  10  jet  diameters  down¬ 
stream.  This  peak  intensi^  was  at  the  same  auqw»ntor  posi¬ 
tion  as  for  the  smaller  nozzle  at  low  exit  velocity  (Figure  12) 

The  possibili^  existed  that  some  of  the  sound  recorded 
was  either  being  reflected  from  or  was  due  to  the  physical 
resonance  of  surrounding  structures.  However,  the  data  pre¬ 
sented  in  Figures  14  and  through  16  indicate  that  these  effects 
were  not  large. 

The  computer  predictions  made  for  the  distributions  of 
axial  veloci^,  turbulence  kinetic  energy  and  pressure  (Figures 
8,  9  and  10)  were  for  the  test  conditions/results  presented  in 
Figure  11.  (RUN  No.  1,  Table  VI.) 

The  peak- intensities  for  the  lower  frequency  sound  noted 
above  occurred  approximately  at  the  predicted  location  within 
the  augmentor  tube  at  which,  (1)  turbulent  mixing  became  essen¬ 
tially  couplets,  (2)  pressure  reached  a  maximum,  and  (3)  the 
axial  velocity  profile  became  fairly  constant.  The  decay  in 
the  predicted  turbulence  kinetic  energy  (Figure  9)  followed 
closely  the  rise  in  noise  amplitude  for  the  two  lower  fre¬ 
quency  ranges  (Figure  11) . 

The  data  were  much  too  limited  to  draw  any  general  con¬ 
clusions.  However,  the  results  indicate  that  the  maximum 
sound  intensities  occur  within  the  augmentor  tube  in  the 
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region  where  jet  mixing  is  nearing  completion.  If  this  ob¬ 
servation  is  correct,  it  would  enable  the  computer  model  to 
be  used  for  full  scale  test  cells  to  predict  the  best  place¬ 
ment  of  sound  suppressors/jet  breaktq>  devices  for  varying 
test  cell  geometries  and  engine  operating  conditions. 


V.  CONCLUSIONS 


Th«  primitiv*  variable  conqputer  model  for  the  turbojet 
test  cell  appears  to  reasonably  predict  the  floiF-  field  for 
engine  eadiaust  conditions  up  to  the  sonic  condition.  The 
conclusion  is  based  upon  subscale  test  cell  msasuremants  of 
velocity  and  pressure  distributions.  To  further  validate  the 
model#  comparisons  should  be  made  with  full  scale  test  cell 
data.  In  addition#  turbulence  Intensity  distributions  mea¬ 
sured  within  the  subscale  test  cell  would  provide  valuable 
new  data  for  additional  model  validation. 
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FIGURE  2.  TURBOJET  TEST  CELL  DIMENSIONS  FOR  NOISE  MEASUREMENTS 
(DIMENSIONS  IN  INCHES) 


FIGURE  3.  AUGMENTOR  PRESSURE  AND  VELOCITY  PROFILES  FOR  STRAIGHT 

PIPE  INLET  AND  ZERO  ENGINE-AUGMENTOR  SPACING,  m  “  0.45  kg/sec 


Me  =  0.69  KgAec 
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\L  DISTANCE  (jet  diometert) 
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FIGURE  10.  PREDICTED  PRESSURE  DISTRIBUTIOMS 


(qp)  BonindiNv 


FIGURE  14.  NOISE  INTENSITY  VS.  FREQUENCY,  RUNS  1  AND  la 
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Lewis  Research  Center 

2100  Brookpark  Road 

Cleveland,  OH  44135 

(Attn:  Mail  Stop  60-6  (R.  Rudley)) 

18.  Federal  Aviation  Administration 

National  Aviation  Facility  Experimsntal  Ctr. 
Atlantic  City,  NJ  08405 

19.  Naval  Air  Propulsion  Center 
Trenton,  NJ  08628 

(Code  PB71:AFK) 

20.  Naval  Ocean  Support  Center 
271  Catalina  Boulevard 
San  Diego,  CA  92152 

(Attn:  M.  Lepor,  M.  Harris;  Code  5121) 
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